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> ■ ABSTRACT 

a^ 

The z = 2.39 radio galaxy 53W002 lies in a cluster of Lya emission line objects 
and may itself be undergoing a major burst of star formation. CO(3-2) emission, at 
0> '. 102 GHz, was detected from 53W002 by Scoville et al. (1997a), who also reported a 

■^ ■ possible 30 kpc extension and velocity gradient suggesting a rotating gaseous disk. In 

Q.. this paper we present new interferometric CO (3-2) observations which confirm the 

previous line detection with improved signal-to-noise ratio, but show no evidence for 
H I source extension or velocity gradient. The compact nature of the CO source and the 

52 I molecular mass found in this object are similar to luminous infrared galaxies and other 

AGNs previously studied. 



1. Introduction 

The distant, compact, narrow-line radio galaxy 53W002, at z = 2.39, has been the subject of 
some interest recently because it may be a relatively young galaxy undergoing a major episode of 
star formation (Windhorst et al 1991), and because it appears to be the dominant galaxy in a field 
of 17 faint Lya emission line objects, of which 7 have been confirmed to be at a redshift close to 
that of 53W002. This is the most distant such cluster known ( Windhorst et al 1991, Pascarelle et 
al 1996a, 1996b), and one motivation for its further study is to test the formation of giant galaxies 
through the assemblage of sub-galactic star-forming clumps (Windhorst et al 1998, Pascarelle et al 
1998). In this respect, analysis of the molecular content and distribution in the dominant galaxy 
of this cluster, 53W002, might bring some new input into the discussion. 

To this end, Yamada et al. (1995) reported a tentative detection of CO(l-O) emission 
from 53W002 using the Nobeyama 45m telescope, which was followed by an interferometric 



measurement of CO(3-2) by Scoville et al. (1997a; hereafter S97a) using the Owens Vahey 
Milhmeter Array. Of particular interest in the latter study was that the CO source appeared to be 
extended by ~ 30 kpc, although the signal-to-noise ratio of the measurement was not very high. In 
other known luminous CO sources, such as Arp 220 and the Cloverleaf quasar (H1413-I-117), the 
molecular gas appears to be concentrated in the inner few hundred parsecs of the galactic nucleus 
rather than being distributed on kpc scales (Scoville et al. 1997b; Kneib et al. 1998; Downes & 
Solomon 1998). As of this writing there are 9 well-documented detections of CO at high redshift 
(recently summarized by Combes et al. 1999), including 53W002, which is one of only two in this 
list that appears not to be gravitationally lensed. 

The presence of abundant molecular gas can be an indicator of widespread star formation 
or nuclear starbursts, and can also be a signpost for the fueling of an active nucleus following a 
merger or interaction between galaxies. Since 53W002 has both an active nucleus and apparent 
extended star formation, as revealed by its extension in the Lya emission line, and therefore 
may be under construction via mergers, observations of molecular lines are clearly germane to 
understanding the system. Given the relevance, and potential uncertainties, of the previous CO 
measurements, we undertook new observations of 53W002 in the CO (3-2) line at 3mm and in 
the continuuum at 3mm and 1.3mm, using the IRAM Plateau de Bure Interferometer, which we 
report here. 



2. Observations and Results 

The observations were performed with the Plateau de Bure Interferometer of the Institut de 
Radio Astronomic Millimetrique over a number of sessions from September 1996 to November 
1997. We searched for the CO(3-2) line emission using two central frequency setups, 101.945 and 
102.195 GHz, for a total of 55 hours. During part of this time ( 20 hours) data at 1.3mm were 
recorded simultaneously to probe the object's continuum at this frequency. The total frequency 
bandwidth was 512 MHz, 50 channels, and typical system temperatures were 130 K at 3mm 
and 500 K at 1.3mm. A number of sources were used to track variations in the phase and gain. 
The source 1739-1-522 was used for flux calibration which could be achieved within an accuracy 
of 5%. The various configurations of the interferometer array, including 4 or 5 15m antennas, 
covered baselines from short (24m) to long (408m). Combining the whole data set results in a 
final synthesized beamsize of 5" x 6" (natural weighting) at PA=90°, and covers a field-of-view 
of 50" . For the results presented hereafter, the velocity origin is at frequency 101.8845 GHz. 
Because of the two different central frequency settings, the (u-v) plane coverage is not identical 
at all frequencies across the CO (3-2) line-profile, which results into a slightly changing spatial 
resolution across the line-profile: for channels centered at -823, -588 and -353 km s^^, the spatial 
resolution is better than for channels centered at -117, -1-117 and -|-353 km s~^. 

CO(3-2) emission was detected with an integrated flux of 1.20 it 0.15 Jy km s^^. The line 
profile is shown in Figure 1 (velocity resolution 118 km s^^), channel maps in Figure 2 (velocity 



resolution 235 km s^^), and the velocity-integrated CO map in Figure 3 (700 km s^^ window). 
The maps in Figure 2 have an RMS uncertainty of 0.45 mJy beam~^, and for Figure 3 the map 
RMS is 0.25 mJy beam~^. The best fit position for the integrated CO source is RA(B1950) = 
17^ 12™ 59^72 and Dec(B1950) = 50° 18' 51.9", with an uncertainty of 0.5" in each coordinate. 
This uncertainty is from two components added in quadrature, one due to phase and baseline 
errors amounting to 0.3", and a second from the map noise of 0.4". For the 235 km s^^ channel 
maps of Figure 2, the positional uncertainty in the channels at —353 km s~^ and —118 km s^^ is 
0.8", and that of the channel at +117 km s^^ is 1.1". 

We find no significant CO emission from any of the Lya galaxies in the field surrounding 
53W002. 

At 3mm (101 GHz) the measured continuum fiux density, using line- free channels, is 0.1 it 0.2 
mJy, yielding a 3a upper limit of < 0.6 mJy. At 1.3mm (230 GHz), using both sidebands, the 
measured flux density is 1.8 it 0.9 mJy. In collaboration with R. Antonucci and T. Hurt, we 
also obtained a 1.3mm measurement using the IRAM 30m telescope 19-channel bolometer of 
1.6 lb 0.5 mJy (February 1997). Combining this with the PdBI result, the weighted average is 
'S'l.Smm = 1-7 ± 0.4 mJy. Any contribution from the centimeter wavelength synchrotron spectrum 
should be negligible at millimeter wavelengths (see Figure 4 of S97a). 

Measured CO and continuum parameters and derived luminosities are summarized in Table 1. 



3. Discussion 

3.1. CO Emission 

First we compare our results with those of S97a. The integrated CO fluxes are in good 
agreement between the two experiments : 1.20 it 0.15 Jy km s~^ found here versus 1.51 it 0.20 
Jy km s~^ found by S97a. Yet, the error quoted by S97a seems to be underestimated for the 
following reason: the RMS given in their Figure 3 is 0.44 mJy beam~^, while the CO line flux is 
integrated over a 5" aperture which comprises two independent beams. Hence the error bar on the 
line flux should read ±0.3 Jy km s^^ (0.44 mJy x ^/2 x 540 km s^^) rather than ±0.2 Jy km s^^ 
as quoted, yielding a 5a detection only. The derived line redshifts between this experiment and 
S97a are consistent within the errors. 

S97a noted a positional shift with velocity in their channel maps, leading to the conclusion 
that the source was extended and possibly rotating. From the velocity-integrated CO map they 
derived a deconvolved CO(3-2) source size of 5.7" x 1.7" (FWHM for an elliptical gaussian), or 
about 30 X 9 kpc (Hq = 75 km s^^kpc^^, qq = 0.5), at position angle 120°. We can neither conflrm 
nor exclude a source of these dimensions, based on our integrated CO map alone. The source can 
equally well be unresolved, 2" in diameter, or 5.7" x 1.7" at any position angle. S97a note that 
"In view of the limited signal-to-noise ratio in the maps, these size estimates are highly uncertain; 



however, it is clear that the emission is resolved, since the centroid is shifted in the channel maps." 

This shift in the centroid is the crucial evidence for extension, but our maps do not show 
it. For direct comparison we plot the same channel maps in our Figure 2 as those of S97a (also 
their Figure 2): —353, —118, +118, and +353 km s^^ (plus we include maps at —588 km s^"*^ 
and —823 km s~^). Between the —118 km s~^ and +118 km s^^ maps of S97a there is a position 
shift of the peak emission of roughly 5", whereas in our channel maps (the three with significant 
emission) there is no shift in the spatial peaks or centroids to within our relative astrometric 
accuracy of ~ 1". Note that emission is also seen at —353 km s"^ in our data, but not in the 
maps of S97a. The lower contours of our integrated CO map (Figure 3) are distorted somewhat 
toward the southwest, but this distortion is not significant within the map noise uncertainties. 
We conclude from our data that a velocity gradient with position such as that seen by S97a is 
not supported, thereby casting serious doubt on the reality of any extension to the CO emission 
in 53W002, at least on the scale of a few arcseconds. The evidence from the Plateau de Bure is 
consistent with a point source spatially coincident at all velocities. 

With regard to the Yamada et al. (1995) CO(l-O) observation, S97a noted that the CO(l-O) 
flux is far too strong to be a counterpart of the observed CO (3-2), given the normal range of CO 
excitation seen in galaxies. The same conclusion applies to our own measurement, which is even 
somewhat weaker. The CO (1-0) as measured by Yamada et al. is an order of magnitude or more 
too strong to be consistent with the CO (3-2) emission. 

The mass of molecular hydrogen in 53W002, as derived from the CO luminosity, depends on 
one's choice of the L'qq -^ M(H2) conversion factor, a. The standard value for Galactic molecular 
clouds based on CO(l-O) data, a = 4Mq (K km s~^ pc~^)~"^, is almost certainly too high for 
dense starbursts and active galactic nuclei. Downes & Solomon (1998) find a value 5 times lower in 
a study of extreme starbursts in ultraluminous infrared galaxies, and Barvainis et al. (1997) derive 
a value 10 times lower in a multi-line study of the Cloverleaf quasar. For L'qq = 1.5 x 10^*^ K km 
s^^ pc^ and using the Downes & Solomon conversion factor, the molecular hydrogen mass would 
be roughly 1.2 x IO^^Mq, while using the Barvainis et al. factor it would be around 6 x IO^Mq. 
The CO luminosities of 53W002 and the Cloverleaf are quite comparable, once account is taken 
for lensing in the latter. If the CO magnification factor is 10 for the Cloverleaf, then its intrinsic 
CO luminosity is L'qq = 1.0 x 10^" K km s^^ pc^ (Barvainis et al. 1997), compared with 1.5 x 10^" 
K km s~i pc2 for 53W002 (Table 1). 



3.2. Astrometric considerations 

In Table 2 we provide a compilation of available astrometric measures of 53W002 at 8.4 GHz 
(Windhorst et al. 1991), 15 GHz (S97a), optical continuum (Windhorst et al. 1991), and CO line 
emission (S97a and current results). The positions summarized in Table 2 show that the location 
of the CO source, which is identical in the three velocity channel maps of the IRAM data set, 



is consistent with that derived from the integrated map in S97a. The CO source location also 
coincides, within the error-bars, with the VLA 8.4 GHz source position (Windhorst et al. 1991), 
while in right ascension it is more than 2" off the position reported for the VLA 15 GHz source 
(S97a). This offset is well beyond the quoted error-bars for the 8.4 and 15 GHz VLA source 
positions. For the current discussion we shall retain the 8.4 GHz positional measure which appears 
to be more precise. A similar offset around 1.7" in right ascension is also found between the CO 
source location and the optical (rest ultraviolet) continuum source (Windhorst et al. 1991). This 
offset however may simply reflect the slight mismatch known to exist between radio and optical 
reference frames and is not considered to be as significant. Indeed, recently Windhorst et al. 
(1998) have noted for 3 other AGNs in the field around 53W002 a mismatch of the same order of 
magnitude between their radio (VLA 8.4 GHz) and optical positions: this reconciles the 8.4 GHz 
and optical continuum positions for 53W002, but, in our view does not explain the disagreement 
between the VLA 8.4 GHz and 15 GHz source locations. 

From this astrometric analysis, we conclude that the CO and 8.4 GHz sources are coincident 
and consistent with the optical continuum source location, according to the above mentioned 
shift between radio and optical reference frames. This leads us to conclude that the CO source is 
related to the AGN located at the heart of 53W002 (Windhorst et al. 1991). 



3.3. Relation of the CO line emission to other components 

As 53W002 has been considered a key object for testing scenarios of galaxy formation, it is 
important to establish its properties solidly. 

An important question for interpreting the evolutionary stage of 53W002 is to determine 
whether the CO line emission is mostly related to the AGN and surrounding activity or to intense 
star formation across the body of the galaxy. Given the fact that the higher signal-to-noise ratio 
IRAM data set in CO does not support the existence of a velocity gradient across the CO line 
emitting region, hence an extension, the validity of associating part of the CO line emission with 
infalling sub-galactic clumps still in the body of a forming galaxy (Windhorst et al. 1998) should 
be questioned. 

On the other hand, taking into account the compact configuration of molecular gas systems 
seen in most ultraluminous infrared galaxies and AGNs, an association of the CO line emission 
with the various components identified by Windhorst et al. (1998) within the central 1" HST 
image seems more plausible. These include the AGN per-se (unresolved source with size ^ 0.06"), 
a compact core with r^ ~ 0.05", a more extended envelope with r^ ~ 0.25" and two "blue clouds" 
roughly colinear across the nucleus, along PA~ 90°, at distance 0.45" W and ~ 0.2" E, which 
are well contained within the 8.4 GHz radio source. The latter also shows an elongation along 
PA= 90° ± 3° on a scale of ~ 0.5". 

Obviously none of the components identified in the HST data can be spatially disentangled 



relative to the CO line contribution. As noted above however, the CO luminosities of 53W002 
and the Cloverleaf quasar are comparable. In the Cloverleaf the molecular gas was found to reside 
in the inner 150 pc region by Kneib et al. (1998), who argued that this is the signature of the 
molecular /dusty torus surrounding the active nucleus. In their early imaging and spectroscopic 
analysis, Windhorst et al. (1991) concluded that the AGN in 53W002 exhibits a continuum 
luminosity close to the boundary with QSOs. By analogy with the Cloverleaf (and ultraluminous 
infrared galaxies in general) it can be argued that the CO line emission detected in 53W002 comes 
mostly from the molecular /dusty torus of the AGN itself. Interferometric data at higher resolution 
are required to resolve this point. 



3.4. Continuum Emission 

Hughes & Dunlop (1998) reported a 3a upper limit of 3.3 mJy for 53W002 in the continuum at 
800 nm. Combining this upper limit (corrected to 800 /im assuming a submillimeter spectral index 
of 3.5) and the value S{CO) = 1.2 Jy km s"^ from Table 1, we find that S'(CO)/58oo > 294 km 
s~^. For four other high-z sources with both CO(3-2) and 800 fiui continuum detections, this 
ratio ranges from 88 km s^^ to 165 km s^^. Therefore, in terms of CO to dust ratio, the case of 
53W002 looks somewhat peculiar. 

Our measurement of 1.7 it 0.4 mJy at 1.3mm is well above the extrapolation of the centimeter 
wavelength spectrum (S97a), suggesting the presence of an extra component. Because upper limits 
only are available for the continuum at 3mm and 850 /um, the nature of this extra component 
cannot be constrained further. If it is dust at about 50K, the upper limit to the mass of dust 
inferred from the measure at 850 /um is Md < 7 x IO^Mq for the cosmological parameters chosen 
in this paper (Hughes and Dunlop 1998). 



4. Conclusions 

This new set of observations of CO(3-2) line emission in the radio galaxy 53W002, collected 
with the IRAM Plateau de Bure interferometer, is of an improved (8 sigma) signal-to-noise ratio 
relative to previous studies. The individual channel maps (averaged over 235 km s^^) do not 
show any positional velocity shift and hence do not support an extension of the CO line emitting 
region beyond 3" (FWHM). Within the quoted uncertainties, the CO (3-2) source is found to be 
positionally coincident with the VLA 8.4 GHz source, as well as with the HST optical continuum 
(rest ultraviolet) image of 53W002. We conclude that a large fraction of the CO line emission is 
likely related to the AGN itself or its immediate surroundings. The CO line profile has a FWHM 
= 420 lb 40 km s~^, while the total CO line flux is found to be 1.20 it 0.15 Jy km s~^. Assuming a 
figure a = 0.4 — O.8M0 (K km s~^ pc~^)~^ for the conversion factor from Lqq to M(H2), which 
better applies to dense regions around the centers of galaxies, we derive a molecular mass in the 



range 6 x 10^ — 1.2 x 10^^ Mq. This value is comparable to the one found within the 150 pc region 
in the Cloverleaf quasar. Therefore, the CO line emission observed in 53W002 could be entirely 
related to the AGN of this galaxy. New measurements of the millimeter continuum of 53W002 at 
1.3mm, as well as upper limits for the continuum at 3mm and 850 /im, imply the presence of an 
extra component in addition to the synchrotron source responsible for the centimeter spectrum. If 
interpreted in terms of dust, then the implied gas to dust ratio in 53W002 is larger than 10^. 

We are grateful to the referee, David Hughes, for very helpful comments, and to the IRAM 
staff for excellent support during the observations performed in service mode. 
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FIGURE CAPTIONS 

Fig. 1. — CO(3-2) line profile, smoothed to 118 km s~^ resolution. Upper scale: velocity relative 
to redshift of z = 2.394 (central frequency 101.8845 GHz). Lower scale: observed frequency. 

Fig. 2. — CO(3-2) channel maps averaged 235 km s~^ velocity ranges. The central cross marks 
the position of the VLA 8.4 GHz radio source (see Table 2). The contours are —3, —2, +2, +3, and 
+4 times 0.45 mJy beam~^ (l'^)- The velocities are relative to z = 2.394, and the 6" x 5" naturally 
weighted beam at PA= 90° is shown. The coordinates are J2000.0. 

Fig. 3. — Velocity-integrated 00(3-2) map of 53W002. Contours range from —3, —2, and +2 
to +6 times the la RMS of 0.25 mJy beam~^. The CO is consistent with a point source (< 3" 
FWHM); the southwest extension of the lower contours is not statistically significant. No emission 
from the Lya emission line objects in the optical field is apparent. The position of the central cross 
is as in Figure 2. The coordinates are J2000.0 



Table 1. CO (3-2) Line and Continuum Parameters 



Parameter 


Value 


Units 


zco 


2.3927 ± 0.0003 




A«FWHM 


420 ± 40 


kni s^^ 


S(CO) 


1.20 ±0.15 


Jy km s^^ 


T r a 


1.48 X 10^'' 


K km s^^pc^ 


S(102 GHz) 


0.1 ±0.2 


mjy 


S{230 GHz) 


1.7 ±0.4 


mjy 



'^ Assumes Ho — 75 km s ^ Mpc ^, 



Qo 



0.5 



Table 2. Positions for various components in 53W002 



Component 


Waveband 


RA(B1950) 


Dec(B1950) 


Reference^ 


CO (IRAM) 


3mm 


17''12"'59=.72±0.05 


50°18'51".9±0.5 


1 


CO (OVRO)'' 


3mm 


17''12'"59=.68 


50°18'51".7 


2 


Radio Continuum 


8.4 GHz 


17''12™59=.76±0.01 


50°18'51".7±0.1 


3 


Radio Continuum 


14.9 GHz 


17''12'"59".86±0.01 


50°18'51".3±0.1 


2 


Optical Continuum 


Gunn g & r 


17''12'"59=.83±0.04 


50°18'51".3±0.4 


3 



''References: (1) This paper; (2) S97a (3) Windhorst et al. (1991) 

''Position estimated from the center of the peak contour of the integrated CO emission from 
Figure 3 of Scoville et al. (1997a). Positional uncertainties unavailable. 
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